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Abstract. Random illumination is proposed to enforce absolute uniqueness and resolve 
all types of ambiguity, trivial or nontrivial, in phase retrieval. Almost sure irreducibility 
is proved for any complex-valued object whose support set has rank > 2. While the new 
irreducibility result can be viewed as a probabilistic version of the classical result by Bruck, 
Sodin and Hayes, it provides a novel perspective and an effective method for phase retrieval. 
In particular, almost sure uniqueness, up to a global phase, is proved for complex-valued 
objects under general two-point conditions. Under a tight sector constraint absolute unique- 
ness is proved to hold with probability exponentially close to unity as the object sparsity 
increases. Under a magnitude constraint with random amplitude illumination, uniqueness 
modulo global phase is proved to hold with probability exponentially close to unity as object 
sparsity increases. For general complex-valued objects without any constraint, almost sure 
uniqueness up to global phase is established with two sets of Fourier magnitude data under 
two independent illuminations. Numerical experiments suggest that random illumination 
essentially alleviates most, if not all, numerical problems commonly associated with the 
standard phasing algorithms. 

1. Introduction 

Phase retrieval is a fundamental problem in many areas of physical sciences such as X-ray 
crystallography, astronomy, electron microscopy, coherent light microscopy, quantum state 
tomography and remote sensing. Because of loss of the phase information a central question 
of phase retrieval is the uniqueness of solution which is the focus of the present work. 

Researchers in phase retrieval, however, have long settled with the notion of relative 
uniqueness (i.e. irreducibility) for generic (i.e. random) objects, without a practical means 
for deciding the reducibility of a given (i.e. deterministic) object, and searched for vari- 
ous ad hoc strategies to circumvent problems with stagnation and error in reconstruction. 
The common problem of stagnation may be due to the possibility of the iterative process 
to approach the object and its twin or shifted image, the support not tight enough or the 
boundary not sharp enough [TH [THl [B5] . Besides the uniqueness issue, phase retrieval is also 
inherently nonconvex and many researchers have believed the lack of convexity in the Fourier 
magnitude constraint to be a main, if not the dominant, source of numerical problems with 
the standard phasing algorithms [U [22j [31]. While there have been dazzling advances in 
applications of phase retrieval in the past decades [25] . we still do not know just how much 
of the error and stagnation problems is attributable to to the lack of uniqueness or convexity. 

We propose here to refocus on the issue of uniqueness as uniqueness is undoubtedly the 
first foundational issue of any inverse problem, including phase retrieval. Specifically we 
will first establish uniqueness in the absolute sense with random illumination under general, 
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Figure 1. Illumination of a partially transparent object (the blue oval) with 
a deterministic (a) or random field A created by a difluser (b) followed by an 
intensity measurement of the diffraction pattern. In the case of wave front 
reconstruction, the random modulator is placed at the exit pupil instead of 
the entrance pupil as in (b). 

physically reasonable object constraints (Figure [T]) and secondly demonstrate that random 
illumination practically alleviates most numerical problems and drastically improves the 
quality of reconstruction. 

To fix the idea, consider the discrete version of the phase retrieval problem: Let n = 
(ni, • • • ,n d ) G 7L d and z = • • • , z<j) G C d . Define the multi-index notation z n = 
z^z^ 2 ■ ■ ■ z n d d ■ Let /(n) be a finite complex-valued function defined on 7L d vanishing out- 
side the finite lattice 

M = jo < n < N| 

for N = (JVi, • • • , Nd) G N d . We use the notation m < n for rrij < rij, Vj. The z-transform 
of a finite sequence /(n) is given by 

F(z) = ^/(n)z- n . 

The Fourier transform can be obtained from the ^-transform as 

F(w) = F{e i27TWl , • • • , e i2 ™ d ) = ^ /(n)e- i27rn ' w , w = (w u ■ ■ ■ , w d ) G [0, l] d 

n 

by some abuse of notation. The discrete phase retrieval problem is to determine /(n) from 
the knowledge of the Fourier magnitude |F(w)|, Vw G [0, l] d . 

The question of uniqueness was partially answered in P, El [13, [18] which says that in 
dimension two or higher and with the exception of a measure zero set of finite sequences 
phase retrieval has a unique solution up to the equivalence class of "trivial associates" (i.e. 
relative uniqueness). These trivial, but omnipresent, ambiguities include constant global 
phase, 

/(•) — > e [d f(-), for some 9 G [0, 2tt], 



spatial shift 

/(•) — > /(• + m), for some m G Z d , 

and conjugate inversion 

/(.)— >r(N-.). 

Conjugate inversion produces the so-called twin image. 

This landmark uniqueness result, however, does not address the following issues. First, a 
given object array, there is no way of deciding a priori the irreducibility of the corresponding 
^-transform and the relative uniqueness of the phasing problem. Secondly, although visually 
no different from the true image the trivial associates (particularly spatial shift and conjugate 
inversion) nevertheless "confuse" the standard numerical iterative processes and cause serious 
stagnation [H [T5J US\> ED] • 

In this paper, we study the notation of absolute uniqueness: if two finite objects / and 
g give rise to the same Fourier magnitude data, then f = g unequivocally. More impor- 
tantly, we present the approach of random (phase or amplitude) illumination to the absolute 
uniqueness of phase retrieval. The idea of random illumination is related to coded-aperture 
imaging whose utility in other imaging contexts than phase retrieval has been established 
experimentally pQ O EJ E, ESS |32] as well as mathematically [28] . 

Our basic tool is an improved version (Theorem [2]) of the irreducibility result of [TTl [T5] 
with, however, a completely different perspective and important practical implications. The 
main difference is that while the classical result [TTl E32] works with generic (thus random) 
objects from a certain ensemble Theorem [2] can deal with a given, deterministic object whose 
support has rank > 2. This improvement is achieved by endowing the probability measure 
on the ensemble of illuminations, which we can manipulate, instead of the space of objects, 
which we can not control, as in the classical setting. 

On the basis of almost sure irreducibility, the mere assumption that the phases or magni- 
tudes of the object at two arbitrary points lie in a countable set enforces uniqueness, up to a 
global phase, in phase retrieval with a single random illumination (Theorem [3]). The absolute 
uniqueness can be enforced then by imposing the positivity constraint (Corollary [T]). For 
objects satisfying a tight sector condition, absolute uniqueness is valid with high probability 
depending on the object sparsity for either phase or amplitude illumination (Theorem [4]). 
For complex- valued objects under a magnitude constraint, uniqueness up to a global phase 
is valid with high probability (Theorem [5J. For general complex- valued objects, almost sure 
uniqueness, up to global phase, is proved for phasing with two independent illuminations 
(Theorem |6| . 

The paper is organized as follows. In Section [2] we discuss various sources of ambiguity. 
In Section [3] we prove the almost sure irreducibility (Theorem [2] and Appendix). In Section 
[2] we derive the uniqueness results (Theorem[3j |4j |5j [6] and Corollary [TJ. We demonstrate 
phasing with random illumination in Section [5] We conclude in SectionjBj 

2. Sources of ambiguity 

As commented before the phase retrieval problem does not have a unique solution. Nev- 
ertheless, the possible solutions are constrained as stated in the following theorem [TTl 127] . 
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Theorem 1. Let the z-transform F(z) of a finite complex-valued sequence {/(n)} be given 
by 

v 

(1) F{z) = az- m Y[F k (z), mGN^eC 

k=l 

where F^, k — 1, ...,p are nontrivial irreducible polynomials. Let G(z) be the z-transform of 
another finite sequence g(n). Suppose \F(w)\ = |G(w)|,Vw G [0, l] d . Then G{z) must have 
the form 

G(z) = |«| e i Vp(n^(z)) [l[F*(l/z*)Y peN^el 

\kei J \fcei c / 
where I is a subset of {1, 2, ...,p}. 

To start, it is convenient to write 

N 

\ f m\ 2 = E E/( m+n )n m K i2 ™' w 

n=-N m&W 
N 

(2) = £ C / (n)e" i2 ™- w 

n=-N 

where 

(3) C/(n)= E/( m + n )/*( m ) 

is the autocorrelation function of /. Note the symmetry Cj(n) = Cf(—n). 

The theorem then follows straightforwardly from the equality between the autocorrelation 
functions of / and g, because F(w)F*(w) = G(w)G*(w), and the unique factorization of 
polynomials (see [27] for more details). 

Remark 1. // the finite array /(n) is known a priori to vanish outside the lattice M , then 
by Shannon's sampling theorem for band-limited functions the sampling domain for w can 
be limited to the finite regular grid 

1 2 2Nj 

2Nj + 1 ' 2Nj + 1 ' " ' ' 2Nj + 1 ' 

since \F(w)\ 2 is band-limited to the set — N < n < N. 



r 1 2 2N- ^ 

(4) M = {(k 1 ,-M--V = i,-,^k 1 = o,-— -,^.} 



There are three sources of ambiguity. First, the linear phase term z~ m in ([!]) remain 
undetermined because the autocorrelation operation destroys information about spatial shift. 
The unspecified constant phase 9 is another source of ambiguity. 

To understand the physical meaning of the operation 

F(z) — ► z- N F*(l/z*) 

consider the case d — 1 

z~ N F*(l/z*) = f*(0)z~ N + f*(l)z l - N + ■■■ + f*(N) 



which is the z-transform of the conjugate space-inversed array {f*(N), f*(N—l), ■ ■ ■ , /*(0)}. 
The same is true in multi-dimensions. 

The subtlest form of ambiguity is caused by partial conjugate inversion on some, but not 
all, factors of a factorable object, with a reducible z-transform, without which the conjugate 
inversion, like spatial shift and global phase, is global in nature and considered "trivial" in 
the literature (even though the twin image may have an opposite orientation). 

In this paper, we consider both types, trivial and nontrivial, of ambiguity, as they both 
can degrade the performance of phasing schemes. Our main purpose is to show by rigorous 
analysis that with random illumination it is possible to eliminate all ambiguities at once. 

3. Irreducibility 
Random illumination amounts to replacing the original object /(n) by 

(5) /(n) = /(n)A(n) 

where A(n), representing the incident field, is a known array of samples of random variables 
(r.v.s). The idea is to first modify the object by the encoding array A(n) so that phase 
retrieval has unique solution and then use the prior knowledge of A to recover /. 

Nearly independent random illumination can be produced by a diffuser placed near the 
object, cf. Figure [TJ The illumination field can be randomly modulated in phase only with 
the use computer generated holograms [5], random phase plates [U [29] and liquid crystal 
phase-only panels [8]. One of the best known amplitude masks is uniformly redundant array 
[12] and its variants [16] . The advantage of phase mask, compared to amplitude mask, is the 
lossless energy transmission of an incident wavefront through the mask. By placing either 
phase or amplitude mask at a distance from the object, one can create an illumination field 
modulated in both amplitude and phase in a way dependent on the distance [52] . 

Let A(n) be continuous r.v.s with respect to the Lebesgue measure on S 1 (the unit circle), 
R or C. The case of S 1 can be facilitated by a random phase modulator with 

(6) A(n) = e^ (n) 

where 0(n) are continuous r.v.s on [0, 2ir] while the case of K can be facilitated by a random 
amplitude modulator. The case of C involves simultaneously both phase and amplitude 
modulations. More generally, A(n) can be any continuous r.v. on a real algebraic variety 
V(n) C C ~ R 2 . For example R and S 1 can be viewed as real projective varieties defined 
by the polynomial equations y = and x 2 + y 2 — 1 = 0, respectively, on the complex plane 
identified as IR 2 . 

The support £ of a polynomial F(z) is the set of exponent vectors in N d with nonzero 
coefficients. The rank of the support set is the dimension of its convex hull. 

Theorem 2. Let {/(n)} be a finite complex-valued array whose support has rank > 2 and 
touches all the coordinate hyperplanes {rij — : j — 1, • • • , d}. Let {A(n)} be continuous 
r.v.s on nonzero real algebraic varieties {V(n)} in C(~ K 2 ) with an absolutely continuous 
joint distribution with respect to the standard product measure on lines "^( n ) where S C N d 
is the support set o/{/(n)}. Then the z-transform of /(n) = /(n)A(n) is irreducible with 
probability one. 
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Remark 2. If the object support does not touch all the coordinate hyperplanes, then the the 
irreducibility holds true, up to some monomial of z. In view of Theorem [7] this is sufficient 
for our purpose. 

Remark 3. The theorem does not hold if the rank-2 condition fails. For example, let p(z) 
be any monomial and consider 

(7) F(z) = £ Cj y(z) 

j 

which is reducible for any Cj G C, except when F is a monomial, by the fundamental theorem 
of algebra (of one variable). Another example is the homogeneous polynomials of a sum 
degree N 



i+j=N 



yj 



which is factorable by, again, the fundamental theorem of algebra. 
The proof of Theorem [2] is given in the Appendix. 

Theorem [2] improves in several aspects on the classical result that the set of the reducible 
polynomials has zero measure in the space of multivariate polynomials with real-valued 
coefficients [TT1 IT8] . The main improvement is that while the classical result works with 
generic (thus random) objects Theorem [2] deals with any deterministic object with minimum 
(and necessary) conditions on its support set. By definition, deterministic objects belong 
to the measure zero set excluded in the classical setting of [T71 [18]. It is both theoretically 
and practically important that Theorem [2] places the probability measure on the ensemble 
of illuminations, which we can manipulate, instead of the space of objects, which we can not 
control. 

In the next section, we go further to show that with additional, but for all practical 
purposes sufficiently general, constraints on the values of the object, we can essentially 
remove all ambiguities with the only possible exception of global phase factor. This decisive 
step distinguishes our method from the standard approach. 

4. Uniqueness 

Without additional a priori knowledge on the object Theorem |2j however, does not pre- 
clude the trivial ambiguities such as global phase, spatial shift and conjugate inversion. For 
example, we can produce another finite array {g(n)} that yields the same measurement data 
by setting 

(9) g(n) = e i V(n + m)A(n + m)/A(n) 
or 

(10) g(n) = e i 7*(N-n + m)A*(N-n + m)/A(n) 

for 6 G [0, 2n] and m G T, d . Expression (JoJ) and (10) are the remaining ambiguities to be 
addressed. 
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4.1. Two-point constraint. One important exception is the case of real-valued objects 
when the illumination is complex- valued (the case of S 1 or C). In this case, on the one hand 
(|9| produces a complex-valued array with probability one unless m = 0,9 = 0,n and, on 



the other hand, (10) is complex- valued with probability one regardless of m. In this case, 
none of the trivial ambiguities can arise. Indeed, a stronger result is true depending on the 
nature of random illumination. 

Theorem 3. Suppose the object support has rank > 2. Suppose either of the following cases 
holds: 

(i) The phases of the object {/(n)} at two points, where f does not vanish, belong to a 
known countable subset of [0, 2ir] and {A(n)} are independent continuous r.v.s on real alge- 
braic varieties in C such that their angles are continuously distributed on [0, 2ir] ( e.g. S 1 or C) 

(ii) The amplitudes of the object {/(n)} at two points, where f does not vanish, belong 
to a known measure zero subset o/R and {A(n)} are independent continuous r.v.s on real 
algebraic varieties in C such that their magnitudes are continuously distributed on (0, oo) 
(e.g. R orC). 

Then f is determined uniquely, up to a global phase, by the Fourier magnitude measure- 
ment on the lattice M. with probability one. 

Remark 4. For the two-point constraint in case (i) to be convex, it is necessary for the 
constraint set to be a singleton, namely the phases of the object at two nonzero points must 
take on a single known value. On the other hand, the amplitude constraint in case (ii) can 
never be convex unless the set is a singleton and the object phases are the same at the two 
points. 



Proof. By Theorem [2] the z-transform of {A(n)/(n)} is irreducible with probability one. We 
prove the theorem case by case. 

Case (i): Suppose the phases of /(ni) and /(n 2 ) belong to the coutable set 6 C [0, 2tt}. Let 
us show the probability that the phase of g(n) as given by Q with m^fl takes on a value 
in 6 at two distinct points is zero. 

Since A(n + m), m^O, and the phases of A(n) are independent, continuous r.v.s on [0, 2ir], 
the phase of g(n),Vn, is continuously distributed on [0, 2tt] for all 9. 

Now suppose the phase of g(no) fc> r some no lies in the set 0. This implies that 9 must 
belong to the countable set 6' which is 9 shifted by the negative phase of /(n + m)A(n + 
m)/A(no). The phase of g(n) at a different location n ^ no, however, almost surely does not 
take on any value in the set B for any fixed 9 e B' unless m = 0. Since a countable union 
of measure-zero sets has zero measure, the probability that the phases of g at two points lie 
in B is zero if m ^ 0. 

Likewise, A*(N — n + m)/A(n), Vm, has a random phase that is continuously distributed 



on [0, 2ir] and by the same argument the probability that the phases of g as given by (10) 
at two points lie in B is zero. 
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Case (ii): Suppose the amplitudes of /(ni) and /(n 2 ) belong to the measure zero set A. 
Since A(n + m),m ^ 0, and A(n) are independent and continuously distributed on R or 
C, the amplitude of g(n) as given by ^ is continuously distributed on M and hence the 
probability that the amplitude of g(n) as given by ^ belongs to A at any n is zero. 



Now consider g(n) given by (10). Suppose that the amplitude of g(n ) belongs to A at 
some n . This is possible only for n = (N + m)/2 in which case g(n ) = e ie /*(n ). The 
amplitude of g(n),n ^ n , has a continuous distribution on R and zero probability to lie in 
A. 

The global phase 9, however, can not be determined uniquely in either case. 

□ 

The global phase factor can be determined uniquely by additional constraint on the values 
of the object. For example, the following result follows immediately from Theorem |3j (i). 



Corollary 1. Suppose that {/(n)} is real and nonnegative and its support has rank > 2. 
Suppose that {A(n)} are independent continuous r.v.s on real algebraic varieties in C such 
that their phases are continuously distributed on [0, 2tt] (e.g. S 1 orC). Then {/(n)} can be 
determined absolutely uniquely with probability one. 

Proof. With a real, positive object, the countable set for phase is the singleton {0} and the 
global phase is uniquely fixed. □ 

4.2. Sector constraint. More generally, we consider the sector constraint that the phases 
of {/(n)} belong to [a, b] C [0, 2ir]. For example, the class of complex- valued objects relevant 
to X-ray diffraction typically have nonnegative real and imaginary parts where the real part 
is the effective number of electrons coherently diffracting photons, and the imaginary part 
represents the attenuation [25J. For such objects, [a, b] = [0,n/2]. 

Generalizing the argument for Theorem [3] we can prove the following. 

Theorem 4. Suppose the object support has rank > 2. Let the finite object {/(n)} satisfy 
the sector constraint that the phases o/{/(n)} belong to [a, b] C [0, 2ir]. Let S be the sparsity 
(the number of nonzero elements) of the object. 

(i) Suppose {A(n)} are independent, identically distributed (i.i.d.) continuous r.v.s on real 
algebraic varieties in C such that their phases {0(n)} are uniformly distributed on [0, 27r] (e.g. 
the random phase illumination ^). Then with probability at least l — \Af\\b — a\^ s ^ 2 ^(2-K)^^ s ^ 2 ^ 
the object f is uniquely determined, up to a global phase, by the Fourier magnitude measure- 
ment. Here [S/2] is the greatest integer at most S/2. 

(ii) Consider the random amplitude illumination with i.i.d. continuous r.v.s {A(n)} C K 
that are equally likely negative or positive, i.e. P{A(n) > 0} = P{A(n) < 0} = 1/2, Vn. 
Suppose \b — a\ < n. Then with probability at least 1 — |A/"|2~^ 5_1 ^ 2 1 the object f is uniquely 
determined, up to a global phase, by the Fourier magnitude measurement. 

In both cases, the global phase is uniquely determined if the sector [a, b] is tight in the sense 
that no proper interval of [a, b] contains all the phases of the object. 

Proof. Case (i): Consider first the expression Q with any m/fl and the [S/2] independently 
distributed r.v.s of g(n) corresponding to [S/2] nonoverlapping pairs of points {n, n + m}. 



The probability for every such the phase of g(n) to lie in the sector [a, b] is \b — a\/(2ir) for 
any 9 and hence the probability for all g(n) with m 7^ 0, 9 7^ 0, to lie in the sector is at 
most I b — o 1 1" 5 / 2 ] (27r) — t 1 ®/ 2 ! . The union over m ^ of these events has probability at most 
\b-a\^(2n)-^ 2 l 



Likewise the probability for all g(n) given by (10) to lie in the first quadrant for any m is 
at most \N\\b-af/ 2 \2Tt)-W 2 \ 

Case (ii): For Q with any m 7^ the [S/2] independently distributed random variables 
g(n) corresponding to [S/2] nonoverlapping pairs of points {n, n + m}, satisfy the sector 
constraint with probability at most 2~^ s ^ 2 ^ if \b — a\ < tt. Hence the probability that all g(n) 
with m/fl satisfy the sector constraint is at most \M\2^ s ^ 2 y 



For (10) with 9 = and any m, g(n ) = f(n ) at n = (N + m)/2 and hence g(n ) lies 
in the first quadrant with probability one. For n 7^ n , g(n) satisfies the sector constraint 
with probability 1/2 if \b — a\ < 71. Now the [(S — l)/2] independently distributed r.v.s 
g(n) corresponding to nonoverlapping pairs of points {n, n + m}, n 7^ n , satisfy the sector 
constraint with probability at most 2~^ s ~ 1 ^ 2 ^ if \b — a\ < it. Hence the probability that all 



g(n) given by (10) with arbitrary m satisfy the sector constraint is at most |A/"|2 1 )/ 2 ^, □ 



4.3. Magnitude constraint. Likewise if the object satisfies a magnitude constraint then 
we can use random amplitude illumination to enforce uniqueness (up to a global phase). 

Theorem 5. Suppose that the object support has rank > 2. Suppose that K pixels of the 
complex-valued object f satisfy the magnitude constraint < a < |/(n)| < b and that 
{A(n)} are i.i.d. continuous r.v.s on real algebraic varieties in C with P{|A(n)/A(n')| > 
b/a or |A(n)/A(n')| < a/b} = 1 — p > for n 7^ n'. Then the object f is determined 
uniquely, up to a global phase, by the Fourier magnitude data on M., with probability at least 

i_I.A/y(*-i)/2]. 

Proof. The proof is similar to that for Theorem |l](ii). 

For Q with any m/fl the [K/2] independently distributed random variables g(n) corre- 
sponding to [K/2] nonoverlapping pairs of points {n, n + m} satisfy < a < |<?(n)| < b with 
probability less than p~^ K l 2 ^ for any 9. Hence the probability that g(n) with m 7^ satisfy 
the magnitude constraint at K or more points is at most \N\p~^ K l 2 \ 



For (10) with any m, |g(n )| = | /(no) I at n = (N + m)/2 and hence g(n ) satisfies the 
magnitude constraint with probability one. For n 7^ no, there is at most probability p for 
g(n) to satisfy the magnitude constraint. By independence, the [(K — l)/2] independently 
distributed r.v.s g(n) corresponding to nonoverlapping pairs of points {n,n + m},n 7^ n , 
satisfy the magnitude constraint with probability at most pK^ -1 )/ 2 !. Hence the probability 



that g(n) given by (10) with arbitrary m satisfy the magnitude constraint at K or more 
points is at most \J\f\p^ K ~ 1 ^ 2 ' 1 ■ 

The global phase factor is clearly undetermined. □ 

As in Theorem [3] case (ii) the magnitude constraint here, however, is not convex. 

4.4. Complex objects without constraint. For general complex-valued objects without 
any constraint, we consider two sets of Fourier magnitude data produced with two indepen- 
dent random illuminations and obtain almost sure uniqueness modulo global phase. 
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Theorem 6. Let {/(n)} be a finite complex-valued array whose support has rank > 2. 
Let {Ai(n)} and {A2(n)} be two independent arrays of r.v.s satisfying the assumptions in 
Theorem [H 

Then with probability one /(n) is uniquely determined, up to a global phase, by the Fourier 
magnitude measurements on M. with two illuminations \\ and A2. 

// the second illumination A 2 is deterministic and results in an irreducible z-transform 
while Ai is random as above, then the same conclusion holds. 

Proof. Let g(n) be another array that vanishes outside M and produces the same data. By 
Theorem [TJ [2] and Remark [T] 

(u \ n(n \ = / ^^/(n + m^A^n + m^/A^n) 

1 1 y{ 1 ' \e w '/*(N-n + m i )\^(N-n + m < )/A < (n), 

for some m; G Z d , 6i G R, i — 1, 2. 

Four scenarios of ambiguity exist but because of the independence of Ai(n),A2(n) none 
can arise. 

First of all, if 

g(n) = e i6i f{n + m^A^n + m^/A^n), i = 1, 2 

then 

e i6l f(n + mi)Ai(n + ra 1 )/X 1 {n) = e i02 /(n + m 2 )A 2 (n + m 2 )/A 2 (n). 

This almost surely can not occur unless mi = m 2 = 0, Q\ — 62 in which case g equals / up 
to a global phase factor. 

The other possibilities can be similarly ruled out: 

g(n) = ^/(n + mOAiCn + mO/AiCn) 

= e ie2 r(N-n + m 2 )A;(N-n + m 2 )/A 2 (n) 

and 

(12) g(n) = e i ^r(N-n + m J )A*(N-n + m i )/A i (n), z = l,2 

for any rrij, 6*,, z = 1, 2. 

The same argument above applies to the case of deterministic A 2 if the resulting z- 
transform is irreducible. 

□ 

5. Numerical examples 

Our previous numerical study (10] and the following numerical examples give a glimpse of 
how the quality and efficiency of reconstruction can be improved by random illumination. 

We test the case of random phase illumination on a real, positive 269 x 269 image consist- 
ing of the original 256 x 256 Cameraman in the middle, surrounded by a black margin (zero 
padding) of 13 pixels in width (Figure [2^a)). We synthesize and sample the Fourier magni- 
tudes at the Nyquist rate (Remark [TJ and implement the standard Error Reduction (ER) 
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(a) (b) 




iidual al each iteration 



20 40 60 80 100 120 140 20 40 60 80 100 120 140 

iteration iteration 

(c) (d) 

Figure 2. ER reconstruction with random phase illumination: (a) recovered 
image (b) difference between the true and recovered images (c) relative change 
\\fk+i — fk\\/\\fk\\ (e) relative residual — |$A/fc| ||/||F|| versus number of 
iterations. 



and Hybrid-Input- Output (HIO) algorithms in the framework of the oversampling method 
[2^1 [25] . By Corollary [TJ absolute uniqueness holds with a random phase illumination. 

Let $ and A be the Fourier transform and the diagonal matrix diag[A(n)] representing 
the illumination. For the uniform illumination A = I. The ER and HIO algorithms are 
described below. 



ER algorithm 

Input: Fourier magnitude data {F(w)}, initial guess Jq. 
Iterations: 

• Update Fourier phase: G k = $A/ fc = \G k (w)\e i&k( - w \ 

• Update Fourier magnitude: G' k (w) = \P-^w)\e' l9k( - w \ 

. Impose object constraint /ttl (n) = { ffg> ^ ™™ * ° 



(a) 



(b) 



e residual al each iteration 



(c) 



(d) 



FIGURE 3. HIO reconstruction with random phase illumination: (a) recovered 
image (b) difference between the true and recovered images (c) relative change 
(e) relative residual versus number of iterations. The final 50 iterations are 
ER. 
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Figure 4. Relative error \\f k 
number of iterations. 
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FIGURE 5. ER reconstruction with uniform illumination: (a) recovered image 
(b) difference between the true and recovered images (c) relative change (d) 
relative residual versus number of iterations. 



In the original version of HIO [13] , the hard thresholding is replaced by 

f (n ^f /fc( n ) iff' k (n)=A-^*G' k (n)>0 
Jk+iW \ f k (n)-Pf' h (n) otherwise 

where the feedback parameter (3 = 0.9 is used in the simulations. ER has the desirable 
property that the residual |||F| — \Gk\\\ is reduced after each iteration under either uniform 
or random phase illumination [131 EH]- When absolute uniqueness holds, a vanishing residual 
then implies a vanishing reconstruction error. 

Figure [2] shows the results of ER reconstruction with random phase illumination. The ER 
iteration converges to the true image after 40 iterations. For HIO reconstruciton we apply 
50 ER iterations after 100 HIO iterations as suggested in [22J . HIO has essentially the same 
performance as ER (Figure 3 (a), (b)). The relative residual curve, Figure [3^d), and the 
relative error curve, Figure 4 however, indicate a small improvement by HIO. The close 
proximity between the vanishing residual curve and the vanishing error curve for ER and 
HIO reflects the absolute uniqueness under random illumination. 
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Figure 6. HIO reconstruction: (a) recovered image (b) difference between 
the true and recovered images (c) relative change (d) relative residual versus 
number of iterations. The final 50 iterations are ER, causing a dip in the 
relative change and residual. 
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Figure 7. Relative error with (a) ER and (b) HIO versus number of iterations. 
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With uniform illumination, ER produces a poor result (Figure [5] (a)), resulting a 72.8% 
error (Figure [5] (b)) after more than 1000 iterations. The relative change curve, Figure [5^c), 
indicates stagnation or convergence to a fixed point after 100 iterations and the relative 
residual plot, Figure ^d), shows non-convergence to the true image. For HIO reconstruction, 
we augment it with 50 ER iterations at the end of 1000 HIO iterations. While HIO improves 
the performance of ER but still leads to a shifted, inverted image which is also severely 
distorted (Figure [6] (a), (b)). The ripples and stripes in Figure [6] (a) are a well known 
artifact of HIO reconstruction [T3J [T5] . As expected, HIO reduces the residual and does not 
stagnate as much as ER (Figure [6] (c), (d)) but its error is greater than that of ER due to the 
interferences from shifted and twin images present under the uniform illumination (Figure 

To summarize, under a random phase illumination, the problems of stagnation and error 
disappear and phasing with ER/HIO achieves accurate, high-quality recovery. These ex- 
periments confirm our belief that a central barrier to stable and accurate phasing by the 
standard methods is the lack of absolute uniqueness. 

6. Conclusions 

In conclusion, we have proposed random illumination to address the uniqueness problem 
of phase retrieval. For general random illumination we have proved almost sure irreducibil- 
ity for any complex- valued object whose support has rank > 2 (Theorem [2]). We have 
proved the almost sure uniqueness, up to a global phase, under the two-point assumption 
(Theorem [3]). The absolute uniqueness is then enforced by the positivity constraint (Corol- 
lary [T|. Under the tight sector constraint, we have proved the absolute uniqueness with 
probability exponentially close to unity as the object sparsity increases (Theorem [4]). Under 
the magnitude constraint, we have proved uniqueness up to a global phase with probabil- 
ity exponentially close to unity (Theorem [5]). For general complex- valued objects without 
any constraint, we have established almost sure uniqueness modulo global phase with two 
independent illuminations (Theorem [6]). 

Numerical experiments reveal that phasing with random illumination drastically reduces 
the reconstruction error, the number of Fourier magnitude data and removes the stagnation 
problem commonly associated with the ER and HIO algorithms. Enforcement of absolute 
uniqueness therefore appears to have a profound effect on the performance of the standard 
phasing algorithms. 

Systematic and detailed study of phasing in the presence of (additive or multiplicative) 
noise with low-resolution random illuminations and sub-Nyquist sampling rates will be pre- 
sented in the forthcoming paper |llj . 

Appendix A. Proof of Theorem [2] 

Our argument is based on [20, 26] and can be extended to the case of more than two 
independent variables. For simplicity of notation, we present the proof for the case of two 
independent variables. 

Proof. First we state an elementary result from algebraic geometry (see, e.g., [33], page 65). 

15 



Proposition 1. If a homogeneous polynomial P(zq, zi, z 2 ) of (total) degree 5 > 2 is irre- 
ducible, then P (zi,Z2) = P(l, zi, Z2) is also irreducible with degree 5. 

For a polynomial Q(zi, z 2 ) of degree 5, the expression 

(13) Q«(z , Zl ,z 2 ) = z 5 Q(-,-) 

defines a homogeneous polynomial of degree 5 with the property Q(zi,z 2 ) = Q^(l, z%, z 2 ). 
The process from Q to Q* is called homogenization while the reverse process is called de- 
homogenization. Homogenization, in conjunction with Proposition [TJ is a useful tool for 
studying the question of irreducibility. 

Let S C N 2 be a given support set satisfying the assumptions of Theorem [2J We now 
show that almost all homogeneous polynomials of 3 variables with the support 

(14) £*• = {{8 - ni - n 2} ni, n 2 ) : n = (m, n 2 ) G S} 
are irreducible. 

We represent the set of all homogeneous polynomials of degree 5 by the projective space 
F Us of dimension i/§ = ( 2+ / ) — 1 • Each homogeneous coordinate of F Us represents a monomial 
of degree 8. The homogeneous polynomials supported on E" are represented by the projective 
subspace 

X={PGP" f : p(n) = 0, Vn ^ E tt } 

where {p(n)} are the coefficients of P. Clearly X is isomorphic to the projective space P s . 
Let Y C X denote the set of reducible homogeneous polynomials supported on We claim 
(cf. [30j, page 47) 

Proposition 2. Y is a closed subset o/X in the Zariski topology. 

A subset of a projective space is closed in the Zariski topology if and only if it is an 
algebraic variety, i.e. the common zero set 

{F 1 = F 2 = ■ ■ ■ = F m = 0} 

of a finite number of homogeneous polynomials F\ , ■ ■ ■ , F m of the homogeneous coordinates of 
the projective space. The Zariski topology is much cruder than the metric topology. Indeed, 
a Zariski closed set is either the whole space or a measure-zero, nowhere-dense closed set in 
the metric topology as stated in the following (see, e.g. [21J, page 115). 

Proposition 3. Any Zariski closed proper subset of a (real or complex) projective variety 
has measure zero with respect to the standard measure on the projective variety. 

Proof of Proposition^ Let the projective spaces P^ J and P^-j represent the homogeneous 
polynomials of degree j and 5—j, respectively, where Uj = ( 2 1" J ) — 1 and us-j = ( 2 t-j J ) — 
Let Yj C Y be the set of points corresponding to polynomials supported on that split 
into factors of degree j and 5 — j. Clearly Y = U^zjYj and we need only prove that each Yj 
is Zariski closed. 

Now the multiplication of two polynomials of degree j and 5 — j determines a regular (i.e. 
polynomial) mapping 
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in the following way. Let G(z) and H(z) be homogeneous polynomials of degrees j and S—j, 
respectively. Let {g(n)} and {h(n)} be the coefficients of G and H, respectively. Then the 
coefficients of the image point H) are given by 

(15) | g{ m ~ n)/i(n) : |m| = m + mj + m 2 = <5 j. 

|n|=5-j 

In other words $ is bilinear in {g(n)} and {/i(n)} and thus is regular. Clearly we have 

Yj = x p^-^ ni 

Since the product of projective spaces is a projective variety and the image of a projective 
variety under a regular mapping is Zariski closed [30J , Yj is a Zariski closed subset of X. 

□ 

Let E = {n 1; n 2 , ■ • • , ns} an d let the ensemble of polynomials corresponding to {A(n)/(n)} 
be identified with 



J] /(n)V(n) = (/(nOVCnO) x (/ (n 2 )V(n 2 )) x ■ • • x (/(n s )V(n s )) C R 2S 

neS 

where 

fV = {{xh - yf 2 , xf 2 + yh) G R 2 : f\ = »(/), / 2 = 3f(/), (x, y) G V}. 

Note that fV is a real algebraic variety in R 2 if V is also. In a similar vein, we now treat 
complex projective space X (resp. variety Y) as real projective space (resp. variety) of 
double dimensions and, by homogenization, we embed Hnes /( n )"^( n ) m ^ an d denote the 
resulting projective variety as V C X. 

Clearly Y fl V is a real algebraic subvariety in V. To show Y fl V is a measure-zero subset 
of V we only need to show that Y fl V C V in view of Proposition [2] and [3j 

Following the suggestion in |20j, we now prove 

Proposition 4. Y n V C V. 

Proof of Proposition [^} It suffices to find one irreducible polynomial in V. 
The argument is based on two observations. First the polynomial 



(16) F(x, y, z) = ax r + by r + cz r 

is irreducible for any positive integer r and any nonzero coefficients a, b, c. This follows from 
the fact that the criticality equations F x = F y = F z = have no solution in P 2 and thus the 
algebraic variety F = is non-singular. 

Secondly, for any S satisfying the assumptions of Theorem |2j there exists a set T C E" of 
three points which can be transformed into {(r, 0,0), (0,r, 0), (0,0, r)}, the support of (16), 
under a rational map. 

We separate the analysis of the second observation into two cases. 

Case 1: (0,0) G E. Then there are at least two other points, say (m, n), (p,q), belonging 
to E. Without loss of generality, we assume p + q = 5. Because E has rank 2, mq — np ^ 0. 
We look for the rational mapping 

(17) z 1 = x kll y k21 z k3 \ z 2 = x kl2 y k22 z k3 \ z = x kl3 y k23 z k33 
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with kij G Z that maps the polynomial 

P{zo, z x ,z 2 ) = c4 + a4- m - n z?z™ + bz\z\ 

to F(x,y,z). This amounts to a linear transformation from the set of independent vectors 

(m,n,5 -m-n), (p,q,0), (0,0,5) 

to the set {(r, 0,0), (0,r, 0), (0,0, r)}. This transformation can be accomplished by the fol- 
lowing matrix 

(m p 0\ r ( qS —n5 —q(6 — m — n) 

n q = — r I — p5 m5 p(5 — m — n) 
S - m -n 5 oymq-npj , Q Q 
mq — np 



where the divisor is nonzero. To ensure integer entries in (18) we set 

r = 5{mq — np) 

and obtain the transformation matrix 

qS —n5 —q(S — m — n) 
—p5 m5 p(S — m — n) 
mq — np 

Case 2: (m, 0), (0, n) G £ for some positive integers m,n. Then there is at least another 
point (p, g)eE such that (m, 0), (0, n), (p, q) are not collinear, which means mn — np — mq ^ 
0. 

Suppose p + q = 5. Consider the polynomial 

P(z , Zl , z 2 ) = az 5 - m z™ + bz^ n z n 2 + cz\z\. 



fell 


ku 


ki3 




h2 


k 2 3 


hi 


k 3 2 


hz 



By the same analysis above the form (16) can be achieved by the transformation matrix 

m p\ r (— q(S — n) q{5 — m) —n(5 — mf 

n q = — I p{5 — n) —p(5 — m) —m{5 — n) 

5 — m 5 — n / [jnn mq np) \ _p n —mq mn 

which has integer entries if r is a multiple of 5(mn — mq — np). With the choice 

r = 5(mn — mq — np) 
the transformation matrix becomes 

— q(6 — n) q(5 — m) —n(5 — m 



hi 


kn 


h 3 


hi 


h2 


k 2 3 


hi 


h2 


hs 



p(S — n) —p(S — m) —m(5 — n) 
—pn —mq mn 

Suppose n = 5. Consider the polynomial 

P(zq, zi, z 2 ) = azl~ m z™ + bz n 2 + czl- v - q z\z\. 



The form ( 16 ) can be achieved by the transformation matrix 

m p \ r I n(n — p — q) q{n — m) n{m — n) 

On q = — I mn — mq — np 

n / nimn — mq — np) \ 

n — mun — p — ql y * \ —pn —mq mn 
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With 

r = n(mn — mq — np) 

the transformation matrix becomes 

f k\\ ki2 k 13 \ I n(n -p-q) q(n - m) 
k%\ k 22 k 23 = mn — mq — np 

K hi h 2 h 3 J \ -pn -mq 

Suppose m = 5. Consider the polynomial 

r[Zo, Zi, Z2) — az 1 + oz Q z 2 + czq z 1 z 2 . 




The form (16) can be achieved by the transformation matrix 



m p \ r I mn — mq — np 

On q = — I pirn — n) mOm — p — q) min — m) 

n I mimn — mq — np) \ 

(J m — n m — p — q I v H \ —pn —mq mn 

With 

r = m{mn — mq — np) 

the transformation matrix becomes 

k\\ k\2 k\z\ ( m n — mq — np 

k 2 \ k 2 2 k 2 3 = I pirn — n) m(m — p — q) m(m — n) 
k 31 k 32 k 33 J \ -pn -mq mn 

To conclude the proof of Proposition |1J in any above case, if the polynomial P(z , Zi, z 2 ) 
is reducible (i.e. has a non-monomial factor), then we can write P = P1P2 and 

(19f (x, y, z) = Pi(z (x, y, z), z x (x, y, z), z 2 (x, y, z))P 2 (z (x, y, z),z x (x, y, z),z 2 (x, y, z)) 

where Pi,P 2 are non-monomial factors. Let / be the lowest (possibly negative) power in 



x,y,z of Pi(zo(x, y } z), z\(x, y, z), z 2 (x, y, z)), i = 1,2. If / > 0, then the factorization (19) 
implies that F{x,y,z) has a non-monomial factor. If I < 0, then the factorization (19) 
implies that [xyz)~ l F{x,y, z) has a non-monomial factor. Either case contradicts the 
that F is irreducible. So P is irreducible. The proof of Proposition [4] is complete. 



act 



□ 



Continuing the proof of Theorem [2| we have from Propositions [3] and [4] that Y fl V is a 
measure-zero subset of V. By dehomogenization and Proposition [T] reducible polynomials of 
a fixed support £ under the assumptions of Theorem [2] comprise a measure-zero subset of 
all polynomials of the same support. The proof of Theorem [2] is now complete. 

□ 
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